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ABSTRACT 



This study employs a radiational model, with large-scale cloud 
parameterization in several layers, in order to compute the absorption 
of solar insolation by the ocean 1 s surface and by atmospheric layers, 
and to compute the planetary albedo for the FNWC primitive equation 
model. Solar insolation dispositions were computed on a gridpoint 
basis using water-vapor mass and cloud amounts. Using emissivity for- 
mulas after Sasamori, the long-wave cooling effects were calculated at 
the earth’s surface and over the tropospheric layers, and were found 
to be dependent upon the cloud parameterization. 

FNWC data over oceanic gridpoints for 16 July 1974 were used to 
test two forms of the cloud parameterization. The objective was to 
determine the parameterization which better verified the radiational 
balance at the top of the atmosphere as a function of latitude as 
compared with satellite climatology. The smaller cloud parameteriza- 
tion gave the better verification. 
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I. INTRODUCTION 



This thesis is a study of a radiative heating parameterization for 
use in the Fleet Numerical Weather Central (FNWC) prediction system. 
The study has as a primary objective, the analysis of the radiational 
and heat balance of the ocean-atmosphere system utilizing FNWC gridded 
data-fields at constant pressure levels for 16 July 1974. The grid- 
point data were selected along four oceanic meridians, three in the 
Pacific and one in the Atlantic, with most of the gridpoints located 
in the Northern Hemisphere (Fig. 1) . 

The specification of amounts of clouds in two designated layers 
has the most influence on the radiative-model dispositions (short- and 
long-wave) . The initial specification of the fractional amounts of 
CL(1) and CL (2) were based on large-scale formulations developed by 
Smagorinsky (1960) and used in a similar study based on the data day 
16 January 1974 by Spaeth (1975) . 

The governing equations of the radiational transfer of the model 
have been derived by Martin (1972, 1974) , who modified both the solar 
and terrestrial radiation transfers to respond to the presence of 
clouds. These clouds were of specified amounts in the two layers, one 
of which is a mid-level cloud and the other a low-level cloud. The 
radiational model in use here has similarities to those in use in UCLA 
and NCAR General Circulation models. 

The physical description of the radiative model may be applied to 
any scale of analysis for which there is adequate resolution of the 
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temperature and moisture data in the vertical. In the horizontal, the 
reliability of the data used here is consistent with that of the FNWC 
analysis to gridpoints, and is typically reported to the nearest tenth 
of a degree with regard to temperature and dewpoint. The radiational 
computations made here are applied to FNWC gridpoints and are designed 
to make a one-hour forward-time step applicable to the FNWC primitive 
equation forecast model, with special adaptations to their G-levels. 

The Smagorinsky (1960) cloud-specification of CL(1) and CL (2) when 
used with the model of this study resulted in global albedo estimates 
that were too high especially in the tropics as compared to the satel- 
lite data of Raschke (1973), thus leading to under-estimates of the net 
radiative balance at the top of the atmosphere and at the ocean surface. 
This could be a result of inadequately resolved relative humidities, 
which in reality might well be lower than inferred by the FNWC tropical 
analysis . 

The comparative results afforded by the 2/3-CL parameterization 
gave reasonably close agreement with the radiative climatology of 
Raschke et al (1973), for the NIMBUS III period 16-31 July 1969. 

Oceanic and tropospheric balances were confuted for both sets of cloud 
conditions (full-CL and 2/3-CL) for 16 July 1974. It was also neces- 
sary to augment the purely radiative model by including a turbulent 
boundary-layer model for sensible and latent heat transports over the 
ocean as adapted from Kaitala (1974) . 

The radiation package results computed with the 2/3-CL cloud model 
compared very favorably with satellite climatology, though only one 
data day (16 July 1974) was used as representative of mid- July sound- 
ings over the ocean, both for the model and the Raschke satellite 
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climatology. However, the latter climatology was based on the data 
period 16-31 July 1969. 

The major effect noted by the use of the 2/3-CL as contrasted with 
full-CL parameterization was the reduction in global albedo and the 
resulting greater contribution to the surface net heating rates. This 
result is in general agreement with recent observations by Von der Haar 
and Hanson (1969) , Von der Haar and Oort (1973) and Raschke et al 
(1973) . However, the major discrepancies in the radiative model re- 
sulted from the specification of layered cloud amounts in the tropics 
(in a manner identical to that used poleward of the tropics) where 
the model gave somewhat excessive albedos. It is felt that the geo- 
metry of solar radiation streams impinging upon subgrid sized cumulus 
cells in the tropics should be remodeled so as to divert a greater per- 
centage of solar radiation downward, as compared with the characteris- 
tic cloud-layering and resultant reflective effects which seem 
permissible in middle and high latitudes. 
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II. DATA PREPARATION 



A. INITIAL DATA FIELDS 

The temperature and humidity data used in this study were arranged 
in the form of soundings taken along four oceanic meridians (Fig. 1) 
of the Fleet Numerical Weather Central (FNWC) Northern Hemisphere ana- 
lyses on 16 July 1974. Computation of radiational dispositions and 
of other heat budget terms were made at these gridpoints. Oceanic 

locations for these computations were chosen because: 

p 

1. The constant a surfaces (where O = — ) of the FNWC primitive 

TT 

equation system are close to being constant pressure levels. 

2. The maritime-area heating rate computations are likely to be 
representative of the month of July 1974 as contrasted with computa- 
tions made for 16 July 1974 for a corresponding set of gridpoints 
over land. 

The three meridians (and their respective number of soundings) 
selected over the Pacific Ocean were located at 125W (25 soundings) , 
170W (25 soundings) and 145E (17 soundings) . The Atlantic Ocean 
meridian was 35W (26 soundings) . This method of selecting "sound- 
ings" along the indicated meridians of the FNWC polar stereographic 
map, made it unnecessary to employ spatial interpolation between 
original data gridpoints along the meridians. Data along line 3 in 
the Pacific was not extended southward of gridpoint (9,55) because 
they fell over land masses (New Guinea and Northern Australia) where 
the surface temperatures and other sounding features were unrepre- 
sentative of the oceanic values. For similar reason two other 
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Figure 1. FNWC polar stereographic grid and meridians (lines 1, 2, 

3, and 4) selected for study. The longitudes A are shown 
for each meridian as well as the extent considered of each 
meridian . 
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gridpoints at (24,24)and (56,6) which also fell just over land masses, 
were replaced with more representative oceanic soundings from the 
gridpoints (23, 24) and (56,7) , respectively. 

The gridpoint soundings in the form of Table I (a) were taken from 
the original FNWC 63-by-63 surface analysis, as well as the vertical 
distribution of T(P) at ten additional levels ranging up to 100 milli- 
bars. Moisture parameters were provided as the vapor pressure at the 
top of the constant flux level, and as the dewpoint depression at five 
selected levels between the 925 millibar and 400 millibar levels. 

These original soundings were then modeled into a radiative sounding 
(Table I (b) ) corresponding in the vertical to FNWC primitive equation 
prediction levels up to 100 mb (see Fig. 2) . 

Data from 0000GMT, 16 July 1974 were used for the Pacific sound- 
ings while 1200GMT, 16 July 1974 data were used for the Atlantic sound- 
ings. These times were used so that the analysis times would correspond 
as closely as possible to local solar noon in each of the indicated 
areas . 

Since the standard instrument level vapor pressure (e ^) was miss- 
ing from the original FNWC analysis at all gridpoints, e . was approxi- 

ai r 

mated using the FNWC field of e^, a computed value of vapor pressure 

at about 20 meters above MSL in the turbulent boundary layer 0 FNWC 

values of T and e over the ocean had been made available using an 
xx 

operational planetary boundary layer model detailed by Kaitala (1974) . 

To perform radiative calculations, it is necessary to have water- 
vapor and C0 2 absorber masses and cloud amounts, CL(1) and CL (2) 
within certain k-level boundaries (Fig. 2). All soundings in this 
study start at sea level with the approximation of surface pressure 
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Table 1(a). Example of typical FNWC sounding for gridpoint (2,2). 
The humidity parameters between 925, ...,400 mb are dewpoint depres- 
sions, whereas the values at the surface and the top of the constant 
flux layer (9999*) are vapor pressures. 



Pressure (mb) 


T(°C) 


Humidity Parameters 


Surface, 1000 mb 


26.1 


e = 23.8 mb 

X 


925 


17.5 


5.4 °C 


850 


13.7 


7.7 


700 


4.8 


12.7 


500 


-11.8 


13.5 


400 


-23.4 


15.3 


300 


-38.2 




250 


-47.1 




200 


-56 o 9 




150 


-67.8 




100 


-81.0 




9999* 


T = 23.5 





x 



* Code 9999 indicates data taken from the top of the constant flux 
level (CFL) of the FNWC initial data program. 
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Table 1(b). Example of the corresponding radiative sounding with mix- 
ing ratio listed at odd k-levels (Fig. 2). Additionally, water-vapor 
and CO ^ absorber masses, and cloud amounts CL(1) and CL (2) (after 
Smagormsky (I960)) are also listed, as these parameters have been 



modeled 


in the radiative 


theory presented 


in this 


study . 


Pressure 

(mb) 


T(°C) 


Mixing 

Ratio 

g/kg 


Absorber Masses 
Water Vapor CC >2 ^ 

(gm/cm 2 ) (cm/cm ) 


Smagorinsky 
Cloud Amounts 
CL (1) , CL (2) 


1000 


26.1 


14.93 








900 


15.7 


8.45 






0.24, CL (2) 


800 


10.9 


0.0 


1.58 


45.53 




700 


4.8 


3 o 04 








600 


-2.4 


0.0 


2.06 


83.53 




500 


-11.8 


1.02 






0.00, CL (1) 


400 


-23.4 


0.36 


2.18 


113.35 




300 


-38.2 


0.15 








200 


-56.9 


0.0 


2.20 


133.99 




100 


-81.0 


0.00 


2.20 


138.67 




0 


-81.0 


0.0 


2.20 


143.35 
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Level- 

Idenlifers 



I 



Layer (0,2) 



U(o), C(0) 
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U(2), C( 2 ) 
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0 = 1.0 



Figure 2. Five- layer radiative sounding used in this study. Levels 
are identified by their values on the k-scale, while layers are identi- 
fied by their level boundary indices in parentheses, e.g. (8,10). 
Pressure-scaled water vapor and CO^ mass increments Au and Ac, respec- 
tively are integrated relative to the surface and the resulting U and C 
are carried at even levels. The temperature T is retained at all levels. 
Amounts of clouds CL(1) and CL (2) in the layers shown have been para- 
meterized for consideration of their radiative effects. 
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7r = 1000 mb. Therefore, the k-levels are associated approximately with 
the FNWC levels P = 1000., 900., 800., 200., 100., 0.0 mb respec- 

ri 

tive to a = 1.0, 0.9, ..., 0.1, 0.0. 

K 

B. INTERPOLATIVE PROCESSING TO K-LEVELS IN RADIATIVE SOUNDINGS 

1. Temperature Profiles 

The gridpoint temperatures were listed at each mandatory level 
of Table 1(a) between 1000., . .., 100 mb. The temperature was assumed to 
be isothermal from 100 mb to p=0.0 mb. The temperature T was set equal 
to the FNWC listed sea-surface temperature. The radiative sounding tem- 
peratures for the remaining k-levels were obtained from either their 
corresponding listed temperature value or by a three-point Lagrangian 
interpolation scheme (Eq. 2-1, (Spaeth, 1975)), to level k when the 
listed temperature-profile does not include the value T^. 

2. Moisture Profiles 

The moisture from the original FNWC soundings (Table 1(a)) were 
converted into mixing ratios at each of the original sounding levels. 

The near-surface vapor pressure (taken as e^) was used to calculate the 
mixing ratio at k=10, as described by Spaeth (1975) : 

q = 621.97 (e /1000) (2-1) 

10 x 

The remaining mixing ratios q were calculated from the original FNWC 
sounding (Table 1(a)), using a dewpoint depression formula (Spaeth, 1975, 
(Eq. 2-3)). The computed q-values were subsequently interpolated to k- 
levels using the previously mentioned three-point Lagrangian procedure 
applied to successive q-values in the original sounding. Resulting 
values of the interpolated q-values are shown for the case of the radia- 
tive sounding at gridpoint (2,2) (Table 1(b)). 
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Due to the fact that most radiosonde humidity data for p<300 mb 
is either unreliable or not available, a power-law extrapolation for- 
mula (Spaeth, 1975) 

V%00 = <V 5 °°> X < 2 - 21 



was used to obtain q-values at k=3,2,l. 

This extrapolative procedure for obtaining q-values at high 
levels was first suggested by Smith (1966), utilizing the least squares 
estimation of X according to 




(2-3) 



with y. = log (q./q c ) and x. = log (p./PrJ* The Smith method was then 
l l b l lb 

further tested by Spaeth (1975) for the vertical scale of the FNWC 

initial data soundings. The resulting correlation coefficient R 

yx 

was found here to lie in the range ,95 to .99, Such high correlations 
indicate that the procedure described by Spaeth (1975) for determining 
upper atmosphere q-values has realistic estimation value for the verti- 
cal scale involved in the radiative sounding 0 
3, Pressure-Scaled Absorber Masses 

The pressure-scaled water vapor absorber mass AU in a layer 
(Fig. 2) were calculated for the five odd k-levels by using the computed 
mixing ratio values (Eq. 2-7, Spaeth, (1975)). The equation for com- 
puting the integrated pressure-scaled water vapor mass is given by the 
algorithm described by Eq, 2-8 of Spaeth. 
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A similar routine was used in the computation of the carbon- 
dioxide scaled absorber mass (Spaeth, 1975, Eqs. 2-10,11). To be noted 
here is the fact that CO^ absorber masses have been stated in terms of 
N.T.P. n pressure-scaled"cms, since the CO^ absorption coefficients are 
generally stated in terms of the reciprocal of this CO^ mass unit. 

C. CLOUD PARAMETERIZATION 

The relative humidities and thus the saturation vapor pressure at 
levels k=5 and k=9 are used in the calculations of the fractional cloud 
cover CL^ and CL^ in layers (4,6) and (8,9) respectively, as depicted 
in Fig. 2. From the computed relative humidity values at k=5 and k=9, 
the fractional cloud amounts for each layer were parameterized using 
the equation (after Smagorinsky, (I960)): 

CL ( 1) = 2.0 (RH ( 5) ) - 0.7 (2-4 (a)) 

CL (2) - 3.33 (RH (9) ) - 2.0 (2-4 (b)) 

Smagorinsky 1 s parameterization of CL(1) and CL (2) permits cloudiness 
fractions of 1.30 and 1.333 respectively with RH=1.0. Values of CL 
greater than 1.0 were considered initially by Smagorinsky (1960) to 
suggest the amount of supersaturation which accompanies precipitation. 

All radiative models in use at the present time limit CL in accord- 
ance with 0 CL < 1.0. In the present radiative package there is little 
opportunity to deduce precipitation and/or supersaturation amounts. 

Hence, the formulations of CL^ and CL^ of Eqs. 2-4 were limited to a 
maximum of 1.0 in each layer in the "full-CL" case. In addition, the 
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(2-5 (a) ) 



Smagorinsky computations for CL^ and CL^ were reduced to 

CL ' = 2/3 (CI^) 

CL 2 f = 2/3 (CL 2 ) (2-5 (b)) 

regardless of whether the original CL values exceeded 1.0 or not, so 
that the revised cloud-cover amounts could be considered for purposes 
of estimating shifts in the global radiative balance. Computations 
made with cloud parameterization of Eq. 2-5, are termed the 2/3-CL case. 

The 1/3 reduction inherent in CL^ 1 , CL 2 ' relative to Smagorinsky 1 s 
CL^, CL 2 were initiated in this study in an attempt to tune cloud 
amounts to give radiation results in closer agreement with recent satel- 
lite climatology of Raschke et al, (1973) . This satellite climatology 
suggests the use of smaller cloud amounts than that specified by 
Smagorinsky. Fractional cloud amounts were considered to be functions 
of large scale effects only. Therefore, small scale convective acti- 
vity, seasonal conditions (except for the soundings) and latitudinal 
effects were not considered in specifying the reduction factor of 
Eq. 2-5. The reduced cloud-parameterization was introduced here for 
estimating large scale radiational effects only. 

D. CLOUD-AREA COVERAGES 

Since Eqs. 2-4a,b or 2-5a,b gave the fractional-coverage of the 
gridpoint area by the appropriate cloud- type, the gridpoint area may 
be thought of as broken into random fractional segments of size 

W (0,0) = (1-CL 1 )*(1-CL 2 ) (2-6 (a) ) 
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wherein there is a combination of clear-over-clear segments in the 
layers. Similarly, the gridpoint area has the fractional area of cloud 
coverage 



W(l,l) = CL * CL 2 (2-6 (b)) 

of upper-cloud amount CL^ overlying lower-cloud amount CL^. Likewise 
the combinations of cloudy-over-clear and clear-over-cloudy, by layers, 
may be visualized as occurring with the weights 

W (1 , 0) = CL 1 * (1-CL 2 ) (2-6 (c) ) 

and 

W (0, 1) = (1 -CL ± ) * CL 2 ( 2-6 (d) ) 

Regardless of whether the full-cloud amounts CL^, CL 2 of Smagorin- 
sky (Eq. 2-4) , or the 2/3-CL amounts of Eq. 2-5 were utilized, it was 
useful to carry, for radiation computations of each sounding, the 
relative weights or fractions of the gridpoint area exposed to the 
specified cloud- layer combinations. Henceforth, the symbols are denoted 
by W ( 0, 0) , W ( 1 , 1 ) , W ( 1 , 0 ) , W (0 , 1) and as given by Eqs. 2-6(a,b,c,d) 
suggest overcast (1) or clear (0) cloud amounts in the indicated layers 
(Fig. 2) , the first index referring to the level of CL^ and the second 
to that of CL 2 * The usefulness of this computational device will be 
exemplified in Sections III and IV, where the procedures for the terres- 
trial and solar radiational computations are discussed and the results 
are summarized over the full set of soundings. 
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III. TERRESTRIAL RADIATION 



A. THEORETICAL AND EMPIRICAL BASIS 

Empirical formulas were developed by Sasamori (1968) for flux emis- 

sivities in the atmosphere as associated with computations required for 

the heat balance requirement of the NCAR General Circulation Model. 

Sasamori derived the empirical emissivity formulas for water vapor and 

CO^ by comparison with the theoretical values built into the Yamamoto 

Radiation Chart (1952) . This chart has proved to be quite accurate for 

computational checking of the Sasamori emissivities and was used by 

Spaeth (1975) and Warner (1974) as a systematic guide for integration 

of the radiative transfer formula developed by Martin (1972, 1974), who 

adapted the Sasamori formulas to the particular layers of interest in 

gridpoint computations of the FNWC primitive equation model (Fig. 2) . 

The essential long-wave flux formulas required for use in the FNWC 

heating package are the following: 

* 

F = net IR flux at earth, k=10 

* 

F _ = net IR flux at level k=6 

6 

* 

F^ = net IR flux at level k=2 

F610 = net IR-flux divergence in the layer (6,10) 

F26 = net IR-flux divergence in the layer (2,6) 

To get the IR-flux divergences in the layers (6,10) and (2,6), the 
* * * * 

differences F - F__ and F - F- must be computed. The detailed 
6 10 2 6 

scheme for making such computations with various combinations of cloud 
cover CL (1 ) and CL (2) is similar to that developed by Martin (1974) and 
as reproduced in detail by Spaeth (1975) . 
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In order to make IR net-flux calculations along the path of integra- 
tion, there must be a physically sound representation of the emissivity 

(G ) as a function of both water-vapor and CC> absorber masses in layers 
wc 2 

along the sounding. For a further discussion of the emissivity formulas 
used in the quadrature scheme, refer to Spaeth's Appendix A (1975). 

B. NET FLUX FORMULATION 
1 . At Level k-10 

The radiative sounding as depicted in Table I (b) , was computed 

as the combination of parameters U(k,10), C(k,10) and T for each re- 

.K 

quired level, k=10 , 8 , . . . , 2 , 1 ,0 . Cloud parameters CL(1) and CL (2) are 
also listed at each gridpoint and in general are non-zero. The grid 
area was then considered to be composed of areal fractions (weights) 
defined in Eqs . 2-6 (a,..., d) and denoted by the symbols W(0,0), W(l,l), 

W ( 1 , 0 ) , W (0 , 1) . 

* 

The composite net flux F^ (CL^, CL 2 ^ at level at each 

gridpoint is then constructed by using the appropriate weight factor to 

★ 

multiply the corresponding reference net flux F computations , defined 

•k k * 

for the special cloud-cover cases, F^ (0,0), F^ Q (1,0), F^ Q (0,1), 

* 

F^q (1,1) . As a result, the composite net flux may be written 
F 10* (CL 1 ,CL 2 ) = W(0 ' 0)F 1 o* (0 ' 0)+W(1 ' 0)F io* (1/0)+W(0/1)F lo : " (0/1) 

+W(1,1)F *(1,1) (3-1) 

Spaeth (1975) has listed these reference net flux formulations in his 
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Eqs • 3-6, 7,8. Using the definitions of W(0,0), W(1,0), W(0,1), W(l,l) 
* 

(CL^ , CL^ ) can be shown to assume the form 

F 10* <CL 1' CL 2> = [l-Ct(2)]{(B 10 - B 6 )-.5[e we (8, 1 0)( B io- B a ) 

+ (e „c (8 ' 10 > +e „c <6 - 10,)( W ]) 

+ (1-CL(2) ) (l-CL(l) ) { b-.5[ (£ (6,10) 

6 wc 

+ e (4,10)) (B -B.) + (e (4,10) (3-2) 

wc 6 4 wc 

+ E (2,10)) (B-Bj + (£ (2,10) 

wc 4 2 wc 

+ G (1,0) ) (B-B )+e ((0,10) ,T ) *B ] } 

wc 2 1 wc 11 

+ CL (2 ) { (B -B ) [1- . 5e (9,10) ] } . 

10 9 wc 

Here 

B, = 1.170403 x 10 -7 T, 4 (3-3) 

k k 



is the Stef an-Boltzmann blackbody flux in ly day at temperature T . 

x. 

Further, e (U , C , 10) is the combined water-vapor and CO 
WC K K 2 . 

emissivity along the path from level 10 to level k. This emissivity is 

considered by Sasamori to be temperature independent for T > 210K, 

whereas represents the temperature dependent emissivity applicable 

for T < 21 OK (see pp. 136-137, Spaeth, 1975). 

* 

The reference net fluxes F^ of Eq. 3-1 are associated with (1) 
clear skies in both layers, (2) overcast in the upper layer only, (3) 
overcast in the lower layer only and (4) overcast in both layers, 
respectively. 

* 



2. Net Flux F 

6 — 

* 

The formula for F (CL ,CL ) has been developed by Spaeth (Eq. 

D 12 

3-10, 1975) in a manner analogous to the derivation of the weighted 
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r io • 
1975) . 



The result is reproduced for completeness after (pp. 48, Spaeth, 



F* = tl-CL(l) ) {B -.5te (6,8) (B — B )+£ (4,6)(B-B) 

6 8 wc 8 6 wc 6 4 

+ (e (4,6)+e (2,6) ) (B -B )+ (E (2,6)+e (1,6 ))(b-b ) 

wc wc 4 2 wc wc 2 1 

+ £ ((0,6) ,T.)*BJ }+(l-CL(l)) (1-CL(2) ) { (B -B ) [1. 

wc 1 1 10 8 

-.5(e (6,8)+e (6,10))]}+CL(1){(B -B,)* 

wc wc 86 

tl-.5e (6,8)]}+CL(l) (1-CL (2) ) { (B -B ) * 
wc 10 8 

[1-.5 (e (6,8)+e (6,10))]} . (3-4) 

wc wc 



* 

3. Net Flux F 2 

The net flux at level k=2 can be calculated in a similar fashion 
★ * 

to F_ and F _ , with the formulation of the individual reference net 
10 6 

* * 
fluxes F^ as previously described. The formula for F^ (CL^CL^) is 

* * 

developed analogous to that for F (CL ,CL ) and F (CL ,CL ) . 

10 1 Z 6 12 

F * = [l-CL(l) ] {B -.5 [e (2,4) (B -B_)+(e (2,4) 

2. 8 wc 4 2 wc 

+ e „c <2 ' 6))(B 6- B 4 )+(E wo <2 ' 6,+e „c (2 ' 8)) * 

(B 8- B 6> +e wc (1 ' 2,(B 2- B l , « TO <<0 ' 2 >' 

+ (l-CL(l)) (1-CL(2))*{(B. -B q ) ( l-.5(e ( 2 , 8 ) 

10 8 wc 

+£ (2,10) ) ] }+CL (1){b-.5[£ (2,4) (B -B ) 

wc 4 wc 4 2 

+ £ (1,2) (B -B-)+e ((0,2), T.)*B.]} (3-5) 

wc 2 1 wc 11 

A new parameter, the total outgoing long-wave radiation to space 

★ 

(FF2) can be readily defined from the expression for F^ by setting to 

zero the terms of Eq. 3-5 representing the downward flux through level 

k=2. These terms are just those denoted by the symbols e ( (0, 2) ,T. ) *B_ 

wc 1 1 
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and € (1,2)* * (B -B.,). If the long-wave radiative model including its 

wc 2 1 

cloud parameterization is realistic, the modeled outgoing flux to space 
(FF2) should compare reasonably with that observed by NIMBUS III satel- 
lite during the same period, after Raschke et al (1973) . The most 
nearly comparable satellite climatological period was 16-31 July 1969. 
These satellite observations were obtained from the NIMBUS III satellite 
atlas of Raschke et al (1973) for the July dates indicated and used for 
comparison with the model values of FF2, point for point (at 5° latitude 
intervals) along the same four meridians of Fig. 1. Also the latitu- 
dinally-distributed mean satellite fluxes were compared with the corres- 
ponding averaged values of FF2. The results showing the comparisons of 
the latitudinally-averaged distribution of the model computations of 
FF2 and the corresponding satellite observations are listed by latitude 
in Table V. 

C. APPLICATIONS TO HEAT BALANCE COMPUTATIONS 

1. At Upper and Lower Boundaries 

In order to compute the heat balance at the top of the earth- 

* 

tropospheric system the composite F^ and the total insolation absorbed 

below level k=2 at each gridpoint are required. Moreover, a radiative 

balance is computed at the surface for each gridpoint using the earth's 

* 

absorbed insolation (Section IV) and the composite net flux loss F . 

In order to consider the possibility of a heat balance at the 
earth's surface, rather than merely a radiative balance, it is necessary 
to add a term representing combined sensible and latent heat transfer 
across the air-sea interface (Section V) . This latter transfer rate at 
each gridpoint was adapted from the FNWC primitive equation model, after 
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the discussion of Kaitala (1974)^ and based upon the FNWC data of 16 
July 1974. 

2. Intermediate Levels 

Consideration of the radiation balance in the atmosphere re- 
quires computation of the long-wave cooling effects caused by the 
flux divergences F26 and F610. F26 and F610 are defined as 

F26 = F 2 * - F 0 * (3-6) 

F610 = F 6 * - F 1q * (3-7) 

where the symbols (2,6) and (6,10) indicate the layer boundaries in- 
volved in the balance consideration for the indicated layer. 

In order to compute F26 and F610 at each gridpoint, complete 

* * * 

listings of IR net fluxes F_ _ , F^ , F_ , associated with each set of 

10 6 2 

reference cloud amounts (0,0), (1,0), (0,1), (1,1) have been computed 

* 

at each gridpoint together with the weighted set of fluxes F^ Q (CL^, 

* * 

CL ), F^ (CL^, CL 2^ ' F 2 < CL i' CL 2 ^ * These cloud-weighted values have 
been constructed using the weighting scheme of Eqs. 2-6, 3-1, 3-4 and 
3-5. A sample gridpoint printout of the net fluxes and flux divergences 
has been included in Table II. 

The determination of the heat balance by layers has been de- 
ferred to Section VI, where the results are displayed in cross- 
sectional form for the 2/3-CL case. 

D. STATISTICAL RESULTS AND COMPARISONS 
1. Net Flux F 1Q * (0,0) 

* 

A statistical test of the F^ (0,0) numerical results was con- 
ducted by utilizing a linear regression program for relating the 
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* 

predictand F (0,0) against the simultaneous values of and 

X 2 = B 10^ as predictors, using all 93 gridpoint data as samples. The 
regression program was based on the BMD02R in the Biomedical set of 
programs (Dixon (1973)). The result was cast in the form of the well- 
known Brunt net flux equation (1932): 

F 10* = B !0 (a+b ^ (3-8) 



where 



B 



10 



St*T 



10 



(3-9) 



is the Stefan-Boltzmann blackbody flux at the surface temperature T^ 
and e is the surface vapor pressure in mb. The multiple correlation 

* i 

coefficient ^l'^2^ WaS 0.99 in the 16 July case, 

as was also true for the data periods of Warner (16 October 1973) and 

Meyers (16 April 1974) , (Table III) . 

* 

Since F may be defined as 



F 



* 

10 



= B 



10 



- F. 



(3-10) 



Eq. 3-10 can then be solved for F n . 

d 

F = B [ (1-a) - bv^] (3-11) 

d 10 

for the clear-sky case. It must be recalled that the computation of 

F. included both the effects of water vapor and of CO^. Since there is 
d 2 

no direct relationship of /e upon the quasi-constant absorber mass of 
CO^, the radiative effect of the CO^ in Eq. 3-11 must be included in 
the term involving (1-a) B^. 
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A cross-seasonal comparison of the coefficients 1-a and -b of 



Eq. 3-11 is presented in Table III, comparing the clear-sky results of 
for 16 July 1975 with those of Meyers and Warner at climatological in- 
tervals of three months earlier and later, respectively. 

The high multiple correlation coefficients indicate that an 
oceanic version of the Brunt downward flux equation is valid for each 
tested season. depends most strongly upon B while only slightly on 

B^v^. This is an apparent result of the nearly constant relative 
humidity over the large oceanic regimes. Moreover, the weak dependence 

of F upon is most evident at the data-time 16 July 1974, when 

d 10 

the coefficient -b is so small as to be statistically insignificant. 

This result is due primarily to the relatively weak gradients of both 

v'e* and in the radiative soundings aloft corresponding to 16 July 1974, 

as compared with the adjacent seasons. 

* 

2. Modification of F ^ for Cloudiness C L^ , CL ^ 

A second statistical test was performed based on a possible rela- 

* 

tionship between the cloudy-sky and the clear-sky cases of F^ . The 

* * 

ratio Y = (CL^ ,CL^ ) (0,0) was used as the predictand in this 

case. The predictor CL was the total opaque cloud cover at each grid- 
point, defined as (after Quinn, 1971 ; and Spaeth, 1975): 

CL = CL (1) + CL (2) (l-CL(l) ) (3-12) 

* 

This regression was run utilizing F (CL^, CL 2 ^ the numerator / com- 
puted using the amounts of CL(1) and CL (2) as calculated from Smagorin- 

* 

sky's formulation and then utilizing F^ (CL^ ' , CL^ ' ) computed from the 
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2/3-CL parameterization. The general form of the regression formula 
required is shown in Eq. 3-13 below, with CL given in Eq. 3-12 

F 10* (CL) = F 10 * ( ° /0) (1 " d(CL)) (3-13) 

A cross-seasonal comparison of the coefficient d for the 2/3-CL and 
full-CL cases using the results of this study and those of Meyers and 
Warner, and of the resulting multiple correlation coefficients are pre- 
sented in Table IV. There seems to be little clear-cut variation in 
the statistics of Table IV between the results of adjacent seasons, 
and little statistical preference between the full-CL and the 2/3-CL 
cases. The high correlation coefficients of Table IV, for all seasons 
considered, indicate the general capability of the IR radiative model 
employed in this study to account for the effective net radiation at 
the surface. 

* 

The mean values of over the 93 gridpoint soundings for 

16 July 1974, as applicable in Tables III and IV were as follows: 
Full-CL Statistics 2/3-CL Statistics 



p * 

10 


(CL) = 


.0774 


ly min ^ 
ly min ^ 


F 10* (CL '’“ 


.0966 


iy 


. -1 
mm 

. -1 
min 


p * 

10 


(0,0)= 


.1449 


F 1q *(0,0)= 


.1449 


iy 


CL = 


.5890 






CL’= .4306 









The superior bar symbol denotes the sample mean. These results indi- 
cate that the surface net flux is decreased by the ratio .4658 with 
CL = .5890, and by the ratio .3333 with CL' = .4306. These results 

indicate that the downward flux F (CL) is 46.6 percent larger at the 

d 

surface than F_(0,0) when full-CL amounts are used in the cloud para- 
d 

meterization model. In the 2/3-CL parameterization the mean downward 
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16 April 1974 


16 July 1974 


16 October 1973 


1-a 


.6557 


.8129 


.6834 


-b 


.0209 


-.0094 


.0161 


F d 


.4527 


.5471 


.4716 


r m 


.9969 


.9950 


.9973 



Table III. Cross-seasonal comparisons of the coefficients 1-a and -b 
from Eq. 3-11, as well as mean values from the Brunt equation and the 
resulting multiple correlation coefficients have been included. 



16 April 1974 16 July 1974 16 October 1973 

2/3-CL full-CL 2/3-CL full-CL 2/3-CL full-CL 



d 


-.73 


-.75 


-.74 


i 

00 

i — 1 


-.72 


-.73 


R 

M 


.9916 


.9945 


.9942 


.9920 


.9974 


.9972 



Table IV. Cross-seasonal comparison of the coefficient d from 
Eq. 3-13 and the resulting multiple correlation coefficients for both 
the 2/3-CL and the full-CL cases. 
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flux at the surface F (CL 1 ) is increased by only 33.3 percent when 

d 

compared to F^(0,0). 

E. COMP ARI SONS OF FF2 WITH SATELLITE CLIMATOLOGY 

Comparison was made of computed-model values of FF2 with satellite 

measurements of total long-wave flux to space (after Raschke et al, 

1973) for the NIMBUS III period of 16-31 July 1969. FF2 is used here 
* 

instead of F^ since the downward flux at k=2 is not measured by the 
satellite. The computation of FF2 was previously discussed. FF2 for 
16 July 1974 was computed for both full-CL and 2/3-CL cases at all 93 
gridpoints. These FF2 values were then interpolated to whole multiples 
of 5° latitude along each of the four meridians and the results were 
then averaged zonally to get a mean latitudinal distribution valid for 
each 5° latitude interval considered in this study. 

Table V shows both of the latitudinally- distributed FF2 cloud- 
model values compared with those extracted from Raschke et al (1973) . 
Raschke' s results were obtained by interpolation from charts using the 
same oceanic meridians as those used in this study. Latitudinally dis- 
tributed total opaque cloud cover CL by Eq. 3-12 are also listed for 
both full-CL and 2/3-CL cases at each latitude. In the bottom line 
of each column in Table V is listed the cosine-weighted mean of each 
set of the column values for the listed latitude range. It should be 
noted that the Raschke results do not correspond to specifically 
known CL values. Table V suggests that the model results for FF2 , 
both full-CL and 2/3-CL, are reasonably close to the corresponding 
values derived from Raschke. A close examination of this table shows 
that in the regions 20S to 5N and 25N to 45N both sets of model results 
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Lat 


RADIATIVE 

FF2 

full-CL 


STATISTICS 

FF2 

2/3-CL 


F2 (RAS) 


CLOUD 
full-CL 
CL Avg 


STATISTICS 
2/3-CL 
CL* Avg 


20S 


.377 


.382 


.40 


.398 


.269 


15 


.386 


.390 


.42 


.456 


.306 


10 


.400 


.403 


.42 


.396 


.258 


5 


.378 


.384 


.39 


.616 


.498 


0 


.357 


.366 


.38 


.615 


.445 


5 


.341 


.354 


.36 


.737 


.534 


10 


.357 


.367 


.34 


.651 


.460 


15 


.369 


.377 


.35 


.622 


.430 


20 


.384 


.389 


.37 


.638 


.431 


25 


.390 


.393 


.40 


.559 


.375 


30 


.372 


.378 


.39 


.601 


.425 


35 


.378 


.382 


.38 


.555 


.383 


40 


.356 


.360 


.36 


.546 


.387 


45 


.333 


.338 


.36 


.671 


.566 


50 


.284 


.298 


.35 


.652 


.475 


55 


.210 


.229 


.33 


.948 


.819 


60 


.269 


.286 


.33 


.635 


.459 


65 


.196 


.219 


.33 


.862 


.590 


Wt 


Avg 


.352 


.360 


.374 


.552 


.405 


Table V. 


Comparison 


of zonally- 


averaged 


longwave flux 


to space, FF2, 



as found by this study for both full-CL and 2/3-CL cloud cases for 
16 July 1974, and by Raschke et al (1973) based upon NIMBUS III measure- 
ments. Also included are composite cloud amount fractions used in the 
two cloud models of this study. Flux values in ly min“l. 
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are slightly less than given by Raschke, while in the region ION to 
20N both sets of values are slightly higher than Raschke. From 20S 
to 45N the agreement between FF2 and F2(RAS) is quite close. From 
50N and higher the correspondence to F2(RAS) is not as good and perhaps 
this is due to the high amounts of cloud cover experienced on this 
day in this region. This could be particularly true since the model 
calculations are for a single day, whereas Raschke ' s calculations 
are based upon 16-day averaged values. Additional limitations of the 
comparisons made here are obvious, when it is recalled that between 
latitudes 20S-5S and between 60N-65N there are fewer than four meri- 
dional lines available for computing the listed zonal values in 
Table V. For all other latitudinal-average values, four meridional 
lines were used in the averaging. While the results from Table V are 
not used to make a final choice of cloud-cover parameterization for 
use in the radiation model, there is evidence of slightly closer agree- 
ment of FF2 with F2(RAS) using the 2/3-CL parameterization. 

■k 

F. COMPARISONS OF CROSS -SEASONAL RESULTS OF F n ^ 

10 

* 

Table VI depicts the latitudinally distributed values of F^ (CL 1 ) 
and of CL 1 at the earth's surface obtained using the 2/3-CL parameter i- 

•k 

zation. A seasonal comparison of the model computed F^ is made with 

data obtained from Meyers (1975) for the 16 April 1974 case and 

from Warner (1974) for 16 October 1973. It is clearly shown that 
* 

F^^ (CL*) is a decreasing function of cloud cover. The most noteworthy 

* 

aspect of Table VI is the belt of high values of F^ for 16 July 1974 
in the latitudinal band 20S to 10S. If the associated CL values are 
examined, it is noted that the composite cloud values were very low 
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16 April 


1974 


16 July 


1974 


16 October 


1973 


Lat. 


F 10* (CL,) 


CL' 


F io* <CL '> 


CL' 


F 10* (CL '> 


CL' 


20S 


.0736 


.842 


.1509 


.269 


.0927 


.628 


15 


.0667 


.838 


.1328 


.306 


.0823 


.629 


10 


.0710 


.752 


.1346 


.258 


.0807 


.607 


5 


.0690 


.647 


.0935 


.498 


.0747 


.707 


0 


.0779 


.659 


.1014 


.445 


.0952 


.489 


5 


.0802 


.591 


.0945 


.534 


.1011 


.471 


10 


.0863 


.639 


.1022 


.460 


.1015 


.459 


15 


.0919 


.548 


.1033 


.430 


.0988 


.486 


20 


.0847 


.536 


.0941 


.431 


.1011 


.437 


25 


.0939 


.473 


.0983 


.375 


.1056 


.383 


30 


.0981 


.446 


.0918 


.425 


.1037 


.448 


35 


.1115 


.520 


.0924 


.383 


.1087 


.484 


40 


.1162 


.408 


.0854 


.387 


.1218 


.325 


45 


.0995 


.467 


.0705 


.566 


.1188 


.301 


50 


.0795 


.674 


.0637 


.475 


.1037 


.381 


55 


.1082 


.555 


.0377 


.819 


.0921 


.452 


60 


.0538 


.816 


.0852 


.459 


.0856 


.450 


65 


.0816 


.796 


.0613 


.590 


.0346 


.885 


Wt 

Avg. 


.0773 


.596 


.0848 


.405 


.0853 


.462 



Table VI. Cross-seasonal comparison of the zonally-averaged net 
fluxes at the surface F and zonally-averaged CL' based on 2/3- 
CL parameterization for the three mid-seasonal dates 16 April 1974, 
16 July 1974, and 16 October 1973. Flux values are in ly min“l. 
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•k 

and this seems to be consistent with the high values for . Further- 

more, the low CL values would appear to be consistent with the northward 

migration of the ITCZ during mid-summer. There is a tendency in each 

* 

season for a relative maximum value of F^ to be located in the sub- 
tropics (latitudes 15-25N) . There is further evidence of a high latitude 

* 

(55-65N) minimum F^ associated with a concentration of maximum cloud 
cover CL*. The analysis does seem to be consistent across seasons. Simi- 
lar conclusions were reached for the full-CL model, although they were 

not included in Table VI. In general, it is felt that the results of 
★ 

the F^q analysis yield no conclusive results in regard to the appropri- 
ateness of the choice of full-CL or the 2/3-CL parameterization models. 
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IV. SOLAR RADIATION 



A. PARTITION OF SOLAR INSOLATION 

The solar constant assumed in this study at level k=0 (top of atmos- 
phere) was 2.00 ly min ^ (Joseph, 1971). This flux was depleted by 4% 
to account for the attenuation caused by oxygen and ozone above the tropo- 
pause. This left the value S = 1.92 ly min 1 at level k=2 to be used 
in this study as the effective solar constant . 

Equation (4-1) was then utilized to compute the effective solar in- 
solation at the tropopause (k=2) as follows: 

F (2) = S [— ] -2 cos z (4-1) 

m 

where S = effective solar constant at k=2 

Cos z = cosine of the zenith angle for the Julian date used 

r/r = ratio of the actual earth-sun distance to the mean earth-sun 
distance for the Julian date used in this study. 

The Smithsonian Meteorological Tables (List, 1958) gives the ratio r/r 

m 

and the solar declination 6 for 16 July 1974, 0000GMT, as these values 
are employed in Eq. 4-1. These values are, respectively: 

— = 1.01644 
r 

m 

6 = 21. 5N latitude 

6 is used in evaluating the cosine of the solar zenith angle, given by 
Cos z = Sin <j) Sin 6 + Cos (p Cos 6 Cos h (4-2) 
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